Transport properties of nano-devices: a one dimensional model study 
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A ID model study of charge transport in nano-devices is made by comparing multi-configuration 
time dependent Hartree-Fock and frozen core calculations. The influence of exchange and Coulomb 
correlation on the tunneling current is determined. We identify the shape of the tunneling barrier 
and the resonance structure of the nano-device as the two dominant parameters determining the 
electron transport. Whereas the barrier shape determines the size of the tunneling current, the 
resonances determine the structure of the current. 
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Nano-systems have been receiving great attention 
lately, mainly because of two reasons. First, they al- 
low the realization of nano-scale electronic devices, such 
as Coulomb blockade structures, transistors, diodes, and 
switching devicesi 1 i 2 i 3 i 4 i 5 i 6 i 7 Second, transport properties 
are governed by quantum mechanics, which opens a novel 
venue for investigating fundamental transport properties 
of quantum few-/many-body systems. The theoretical 
description of quantum transport properties has proven 
to be a true challenge. Conventional theories, such as 
the Landaucr-Buttiker approach and density functional 
theory do not perform wcll, 8 - 9 - 1( iii*i£ and are orders of 
magnitude off the experimentally measured currents. Re- 
cently, configuration-interaction (CI) analysis was shown 
to give a much better agreement^ with currents about 
three times the experimental values. 

Despite this progress, the physics underlying electron 
transport in nano-devices is not yet well understood. 
Here, a first step is made to close this gap. Our analysis 
is based on a one dimensional multi-configuration time- 
dependent Hartree-Fock (MCTDHF) method^. The 
MCTDHF approach takes full account of exchange and 
correlation effects. Whereas in CI the basis is fixed and 
only the expansion coefficients are optimized, in MCT- 
DHF both coefficients and basis functions are optimized. 
As a result, a much smaller basis set can be used. This 
proves important especially for transport processes in 
continuum- like structures, where a CI analysis requires 
a large number of basis functions for convergence. 

The influence of many-body effects on transport phe- 
nomena presents a central issue in the characterization of 
nano-devices. We use the following approach to investi- 
gate this question. The MCTDHF analysis is compared 
to time-dependent Hartree-Fock (TDHF), frozen-core 
Hartree-Fock (FCHF) and frozen core Hartree (FCH) 
calculations. Comparison of the different approaches re- 
veals the role of exchange and correlation in the electron 
transport. It also gives a sense of the quality of various 
approximations used in the electron transport analysis 
of nano-devices. Further, a correlation measure is intro- 
duced with the goal to quantify the influence of many- 
body effects. Finally, we identify two parameters that de- 
termine the tunneling characteristics, namely, the shape 
of the barrier and the resonance structure of the nano- 



device. 

A schematic of the ID model potential U(x) of the 
nano-device is given in Fig. ^ The potential is deter- 
mined by U(x) = Uq for |x| < I and U(x) — V n d{x) + 
V c (x) otherwise. Atomic units are used throughout the 
paper, where 1 bohr — 0.0529 nm, 1 hartree = 27.2 eV, 
and io — 6.6 x 10 -3 A is the atomic unit of current. The 
parameter Uq represents the depth of the nano-device and 
its length is 21. We assume that the device is neutral, i.e., 
it contains the same number of ions and electrons, f c = 4. 
Therefore, for |a;| > I, U(x) is composed of the ion po- 
tential V n d{x) = — YjiLi V\/ a n + ( x ~ di) 2 , and of the 
effective conductor potential, V c (x) = Uf[2 — tanh((d + 
x)/w) — tanh((d— x)/w)]/2. Here di is the position of the 
ith nucleus, a n is the shielding parameter of the ID ion 
potential, and Uf is the effective potential in the conduc- 
tor. Parameters d and w are used to model the potential 
in the transition region between nano-device and metal 1 -. 

Further the Coulomb interaction between the / = 
/c + 1 electrons is modelled by (1/2) Y^,{j Vee(%i>2j) with 
V ee (x i: Xj) = 1/ ^Ja 2 + (xi — Xj) 2 . The shielding param- 
eter a removes the Coulomb singularity. As the averaging 
over the dimensions transversal to the direction of elec- 
tron transport results in an effective, shielded Coulomb 
interaction, a can be viewed to represent the transversal 
properties of the nano-devices*. 

For all the calculations presented here the number 
of electrons is fixed to / = 5, four core electrons in 
the nano-device plus one electron tunneling through the 
device. The wavefunction of the transport electron is 
given by a Gaussian initial wavepacket, ip t , with 1/e- 
width a, momentum po, and energy Eq — Pq/2, see 
Fig. [IJ. The initial distance from the tunneling device, 
xo = —200 bohr, is chosen large enough so that the total 
multi-electron initial wavefunction can be written as an 
anti-symmetrized tensor product of the electronic ground 
state of the nano-device, f c , and of the transport elec- 
tron, i.e., * = A[^ c (qi,q 2 ,q3,qA) ® Vtiib)]- Here, A de- 
notes the anti-symmetrization operator, and = (xi, Si) 
{i = 1, 5) the space-spin coordinate. 

The MCTDHF current is a combination of various sin- 
gle and multi-electron processes. To untangle the differ- 
ent contributions and their importance we calculate the 
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FIG. 1: The effective ion potential U(x) (thick solid), the ini- 
tial Gaussian wavepackect (thin solid) , and the effective frozen 
Hartree potential u F f f H (x) = U(x) + u FCH (x) (dashed). The 
Coulomb potential u FCH (x) is the Hartree contribution of the 
initial core. 



electron transport by different methods. In the MCT- 
DHF approach the exact wavefunction is approximated 
by the ansatz 



* = -™ E A h-j f t)-Vj, (<?/,*)■ (1) 

Note that the number of wavefunctions n in the ansatz 
is larger than the number of electrons, resulting in 
Slater determinants. The ansatz is inserted into the 
Schrodinger equation and the Dirac-Frenkel variational 
method is used to derive a coupled system of time- 
dependent, nonlinear, partial differential equations for 
the coefficients A and wavefunctions tp. For more in- 
formation on the MCTDHF method, see Ref. Q In 
the calculations the simulation interval is chosen from 
x = —400 bohr to x = 400 bohr with the number of grid 
points N — 800. Convergence is obtained for n = 11 
wavefunctions. Increasing n and N changes the electron 
density by less than 0.1%. 

The second method used in our investigation is time- 
dependent Hartree Fock (TDHF), which is MCTDHF 
for f — n. The third and fourth methods are frozen 
core methods, where the correlated ground state wave- 
function of the nano-device is calculated by MCTDHF 
and assumed to remain frozen during the transport of 
the conducting electron. In dependence on whether ex- 
change effects between ground state and conducting elec- 
tron are taken into account, these approaches are termed 
frozen core Hartree (FCH) and frozen core Hartree Fock 
(FCHF). The frozen core approximation reduces the 
multi-electron Schrodinger equation to a single-electron 



equation given by 

idt<f>(x, t) = h<t>(x, t) + J V ee [x' \x)p c {x')dx' <j>{x,t) 
- J rW{q'\q)V ee (x', x)(j>{x', t)6 s , s ,dq' 
-J rM(q'\q)h(x')cl>(x',t)5 s , s ,dq' 

-2 J T^(q',q"\q , ,q)V ee (x',x")<l>(x",t)S S!S , l dq'dq", 

(2) 

where s is the spin index of the conduct- 
ing electron, Pc{x) = J Tc (q\q)ds, and 
h = p 2 /2 + U(x). The core electron reduced 

density matrices are defined as if^o^ |gi) = 
[2/(/ c -l)] / lf 3 ( ? i, q 2 \ qi , q 2 )dq 2 and if q' 2 \ Ql , q 2 ) = 

( 2 °) J **0?i, q' 2 , 93, g/J*c(?i, «2, 93, ...,q fc )dq 3 ...dqf c 17 . 
The first two terms on the right hand side of Eq. J2J give 
the frozen core Hartree (FCH) approximation. The last 
three terms account for the ground state correlation of 
the nano-device and the exchange effect between bound 
and tunneling electrons. In order to keep the problem 
as simple as possible, we focus here on the field-free 
scattering of the electron from the nano-device. In the 
weak field limit, where the field induced distortion of the 
bound state is weak, the change of the initial energy of 
the incoming electron Eq is equivalent to applying a bias 
voltage V ~ Eg determined by the chemical potential 
of the two leads. One atomic units of energy change 
corresponds to a bias change of 27.2V. 

In Figs. |2 and 13 we have plotted the average current / 
flowing through the nano-device as a function of the en- 
ergy of the incoming electron. The parameters have been 
chosen to reflect experimental systems, see the figure cap- 
tion. The device size of around 0.5 nm, the ionization po- 
tential of 10 eV, and the Fermi energy of — 5eV are typi- 
cal values found in experimental nano-molecule systems, 
such as in the Au/dithiolated-benzene/Au junctionAi. 
The average current I(x) =< i(x,t) > is calculated be- 
tween the initial time and the time when the reflected 
pulse returns to its initial position. The time dependent 
tunneling current i(x,t) is obtained by 

i(x,t) = — 
ii 

n 

x^2pj. k [f*{x,t)d x (p k {x,t) - <f* k (x,t)d x (pj(x,t)], (3) 
j.fe 

where p j>k = YlZ,-if A j,i 2 -i f A k,i 2 -if denotes the den- 
sity matrix and I(x) is calculated at x ~ 80 bohr. The 
electron energy Eq is measured with reference to Uf. 

The major difference between Figs.|21and|3|is the use of 
different shielding parameters a — 4.5 bohr and a = 3.0 
bohr, respectively. The decreasing value of a reflects 
an increasing electron-electron interaction strength. Al- 
though it is difficult to assign the ID Coulomb correlation 
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strength resulting from different values of a to specific de- 
vices, some qualitative considerations can be made. We 
believe that a ss 2 to 5 represents best the 3D experimen- 
tal situation. On the one hand, calculations for a = 7 
show hardly any interaction between the transport and 
the core electrons. On the other hand, for a < 2, the 
electron repulsion becomes too strong. As a result, the 
ionization potential increases and can no longer be fitted 
to values typical of experimental devices. 

The MCTDHF, TDHF, and FCHF results are depicted 
by the full, dotted, and dashed lines, respectively. The 
dash-dotted line in Fig. [5] represents the FCH result. 
All curves show the same general trend. Whereas the 
structure of the curves is determined by the resonances, 
their growth with increasing Eq depends mainly on the 
height and thickness of the tunneling barrier. There are 
two resonances in the ranges Eq — 0.02 ~ 0.04 and 
Eq = 0.1 ~ 0.14. The strength and width of the second 
resonance is more pronounced mainly due to the weaker 
tunneling barrier for higher Eq. 

The comparison of the various approaches reveals the 
role of exchange and correlation effects in determining 
the current through the model nano-system. The FCH 
approach neglects two important effects: (i) exchange in- 
teraction between conducting and bound electrons; (ii) 
interaction between tunneling and core electrons in the 
sense that the tunneling electron induces a perturbation 
to the bound electrons of the nano-system. The result- 
ing potential change modifies the tunneling current. In 
order to capture this effect, the frozen core approxima- 
tion has to be removed. In the FCHF approximation the 
exchange interaction is taken into account, whereas ap- 
proximation (ii) remains. In the TDHF approach part of 
the limitation (ii) is removed by making the wavefunc- 
tion, that governs the bound electron dynamics, time- 
dependent. However, TDHF does still not fully account 
for the cross-talk between conduction and bound elec- 
trons, as correlation is neglected by the product ansatz 
of single-electron wavefunctions. Finally, in MCTDHF 
all effects are treated properly. 

Figure [21 shows that FCH introduces a significant error 
with regard to the position of the resonances, as can be 
seen from the first resonance. As a result of the large 
width of the second resonance, the difference between 
FCH and MCTDHF is blurred. In Fig. FCH was 
not plotted as the error in resonance position becomes of 
the order of the difference between first and second reso- 
nances. Comparison of FCH and FCHF in Fig. [3 shows 
that a substantial part of the resonance shift comes from 
the neglect of the exchange effect. 

The remaining difference can be attributed to the per- 
turbation of the electrons in the nano-device caused by 
the tunneling electron. The resulting change of the elec- 
tronic core modifies the position of the resonances and 
therewith the tunneling current. A comparison of the 
FCHF, TDHF and MCTDHF in Fig. demonstrates 
the importance of correlation in the interaction between 
bound and tunneling electrons^. Although the TDHF 



0.001 




0.02 0.04 0.06 0.08 0.1 0.12 0.14 
E (a.u.) 

FIG. 2: Average current versus energy Eo, obtained by MCT- 
DHF (full), TDHF (dotted), FCHF (dashed), and FCH (dash- 
dotted). Parameters are: f = 4, a = 4.5 bohr, 21 — 10 bohr, 
d = 11.0 bohr, w = 0.8 bohr, a = 40 bohr, x = -200 
bohr, U = -0.92 hartree, U f = -0.18 hartree (-5 eV), and 
I p = 0.382 hartree (10.4 eV). 
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FIG. 3: Average current versus energy Eo, obtained by MCT- 
DHF (full), TDHF (dotted), and FCHF (dashed). Parameters 
are the same as those given in Fig. [5] except for a = 3.0 bohr. 
The following parameters had also to be modified to keep the 
ionization potential I p — 0.388 hartree (10.6 eV) and reso- 
nance structure as close as possible to the system in Fig. 2: 
21 = 12 bohr, and U = -0.97 hartree. 



approach performs better than FCHF, it is in most places 
off by an order of magnitude from the MCTDHF calcu- 
lation. This highlights the fact that during tunneling a 
highly correlated state between core and transport elec- 
trons builds up. The amount of correlation strongly de- 
pends on the electron-electron interaction strength, as 
can be seen from a comparison between Figs. and |21 

Finally, over most of the energy range (Eq < 0.12) the 
electron transport is elastic. The nano-device returns to 
its ground state after the conducting electron has passed. 
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FIG. 4: Dashed lines: correlations of the initial bound states, 
n c (a,t — 0). Solid curves: overall correlations measured at 
the tunneling peak of the nano-devices, n? k (a), the peak value 
of n c (a,t). The difference between solid and dashed lines is 
the contibution due to the interaction between tunneling and 
core electrons. 



Only for energies (E > 0.12) we find evidence of excited 
state population. The excitation causes a suppression of 
the resonance at E = 0.12 in Fig- El This " strong- field" 



limit is not analyzed further here. 

One way to characterize the build up of Coulomb 
correlation between tunneling and bound electrons is 
to consider the occupation number of the natural or- 
bitals, which are the eigenvalues of the density matrix 
/>it The number of natural orbitals is n > f and they 
are ordered with respect to their eigenvalue. The low- 
est orbital refers here to the largest occupation num- 
ber. In the uncorrelated HF limit only the lowest / or- 
bitals are occupied. Therefore the Coulomb correlation 
can be measured by n c (a,t) = ^i( a t~t)/ ft which 

is the occupation number of the natural orbitals with 
i > f. Here, S"=i A >(M) = / MM) > A i+ i(a,t), and 
0<n c (a,t)< 1-f/n. 

The correlation measure is plotted in Fig. 0] for the 
parameters of Figs. [21 and [31 Up to Eq w 0.06 the cor- 
relation introduced by the tunneling electron is smaller 
than the initial core correlation n c (a, t = 0), which is in- 
dependent of the value of E . In this range the physics 
is dominated by the bound electrons that determine the 
position of the resonances. For larger energies, the inter- 
action between bound electrons and tunneling electron 
adds a significant amount of correlation to the system. 
The sharp increase of n c (a,t) agrees with the large dif- 
ference between MCTDHF and the other approaches in 
Figs. Hand El 
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